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ABSTRACT 

The dynamical fluctuations inside the quark and antiquark jets are stud- 
ied using Monte Carlo method. Quark and antiquark jets are identified from 
the 2-jet events in e + e~ collisions at 91.2 GeV by checking them at parton 
level. It is found that transition point exists inside both of these two kinds of 
jets. At this point the jets are circular in the transverse plane with respect to 
the property of dynamical fluctuations. The results are consistent with the fact 
that the third jet (gluon jet) was historically first discovered in e + e~ collisions 
in the energy region 17-30 GeV. 
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1 Introduction 



As is well known, in high energy collisions the hadronization of partons, being a soft 
process, can not be analyzed using perturbative QCD (pQCD). Therefore, phenomeno- 
logical analysis has to be applied to study the properties of the final state hadronic 
system, in order to get information which would lead to a better understanding of the 
dynamics of strong interaction. Investigating the dynamical fluctuations inside jets 
produced by the fragmentation of partons is an effective way to explore the dynamics 
in high energy collisions. 

In a 2-jet event in e + e~ collisions, the two produced jets are expected to be back 
to back and the thrust (or sphericity) axis coincides with the jet axes. Furthermore, 
the two jets being produced from a quark and an antiquark respectively are expected 
to have similar property. Under these assumptions the study of the dynamical fluctua- 
tions inside a single jet could be performed through the investigation of the fluctuations 
in a 2-jet event. Recently, this kind of investigation has been carried out both theoret- 
ically IJ and experimentally 0. 

In view of the importance of this problem, it is worthwhile studying the dynamical 
fluctuations inside a single jet directly |J. The aim of this paper is to carry on this 
study for quark and antiquark jets separately using Monte Carlo method. Through this 
study, the influence of the above mentioned two assumptions are gotten rid of and the 
dynamical fluctuation properties in the hadronization process of quark and antiquark 
can be directly obtained. 



2 Selection of quark and antiquark jets 

Event samples of e + e~ collision at cm. energy y/s = 91.2 GeV are produced using 
LUND JETSET7.4 Monte Carlo code. Only charged particles are used in the calcula- 
tion. 

Two-jet events are selected using the Durham jet algorithmic . In this scheme, a 
jet resolution variable y cut is defined for every pair of particles (or jets) i and j in an 
event by: 

2mm(Ef,E]) H 
Vij = — (1 - cos%) (1) 

where Ei and Ej are the energies of the two particles (or jets), % is the angle between 
them and s is the square of cm. energy of the event. Jets or particles with < y cut 
are combined into a single jet. This procedure is repeated untill all pairs i and j satisfy 
Vij > Vcut- The 2-jet events remained at the end of this process consist the 2-jet sample 
corresponding to the cut parameter y cut . 

Since we are using Monte Carlo method, it is possible to study the 2-jet events both 
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at parton level and at hadron level. Quark jet and antiquark jet can then be identified 
through matching the jets in these two levels using the following technique H]. 

Firstly, identify the type of a single jet at parton level by checking every parton 
involved in this jet and giving a weight +1 to quark, -1 to antiquark and to gluon. 
The sum of the weights is expected to yield +1 for quark jet and -1 for antiquark jet. 
The cases where the sum of weights in a jet is neither +1 nor -1, or the weights in the 
two jets of one event are both +1 or both -1 are neglected. This cut condition throw 
away about 1% of the total events. 

Secondly, consider the configuration of the 2-jet event at parton and hadron levels. 
The directions of hadron jets might not be exactly that of the corresponding parton 
ones, since the effects of hadronization may change their directions. 

The type of a hadron jet can be determined by checking the parton jet closest to 
it. The hadron jet closest to the parton level quark jet is considered as the quark jet 
and that closest to the parton level antiquark jet is identified as the antiquark jet. In 
this way, the quark jet subsample and the antiquark jet subsample are obtained. 



3 Single-jet coordinate system 

In the 2-jet events, the thrust (or sphericity) axis is supposed to be the direction in 
which the quark- antiquark move back to back with high momenta. Therefore, it is 
chosen as the longitudinal direction in the physics analysis. While when the single 
quark or antiquark jet is considered separately, the thrust axis is no longer appropriate 
to be considered as the longitudinal direction. The direction of the total momentum of 
the particles inside the jet should be taken as the longitudinal axis instead. Thus a jet 
coordinate system with the third axis pointing along the direction of jet momentum 
should be built up. The transformation from the lab system to the jet system is 
acomplished through the following two steps. First, rotate the lab frame o-xyz an 
angle (p around the z axis to form the o-x'y'z' frame and then turn the latter an angle 
9 around y 1 axis to get the single-jet coordinate system o-x"y"z", cf. Fig.l. This 
transformation can be expressed as: 

(x" \ I cos 9 cos (f cos 9 sin ip — sin 9 \ / x \ 
V" = -sin^ cos^ L (2) 

z" / \ sin 9 cos ip sin 9 sin ip cos 9 J \zj 
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Fig.l The transition from laboratory frame to quark (antiquark) jet frame 
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Fig. 2 The distribution of the angle between the quark (antiquark) jet axis 
and the thrust axis, (a) quark jet (b) antiquark jet 

The distribution of the angular difference between the direction of the quark (an- 
tiquark) jet axis and the thrust axis is shown in Fig.2. It can be seen from the figure 
that the jet axis is very close to the thrust axis. 



4 Dynamical fluctuations 

The dynamical fluctuations are characterized by the anomalous scaling of normalized 
factorial moments (NFM) |6j 

FgjM) = — J2 {nminm ~ ' >m " q + - ^ M^oo, (3) 

where a region A in one-, two-, or three-dimensional phase space is divided into M 
cells, n m is the multiplicity in the mth cell, and (■ ■ •) indicates vertically averaging over 
the event sample, consisting of quark or antiquark jets. 

As Ochs pointed out, the dynamical fluctuations occured in higher- dimentional (2D 
or 3D) phase space have the projection effect on the fluctuations in lower-dimentional 
space causing the second-order ID NFM goes to saturation by the rule f|: 

F 2 (a) (M a ) =A a - B a M~ 7a , a = y,p t ,<p (4) 
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where the exponent 7 a describes the rate of going to saturation of the NFM in the 
direction a and is the most important characteristic for the higher-dimensional dynam- 
ical fluctuations. If the values of 7 in two directions are equal, j a = 75, the fluctuations 
are isotropic in the a, b plane; while when 7 a 7^ 7& the fluctuations are anisotropic in 
this plane || . The nature of the dynamical fluctuations (or that of the fractal) can be 
expressed in terms of the Hurst exponent H a i,, which can be obtained from the values 
of 7 a and 75 as [g 

H ab = — 5 

1 + la 

The dynamical fluctuations are isotropic (self-similar fractal) when H ab = 1, and 
anisotropic (self-affine fractal) when H ab / 1 0. 

In our Monte Carlo simulation a total number of 2000000 events of e + e~ collisions at 
91.2 GeV are produced by JETSET 7.4 generator. Different y CVLt are used to select the 
2-jet events in the phase spase region [0 < y < 5], [0.1 < p t < 3.0], [0 < <p < 2n] and the 
second order ID factorial moments F 2 (y), F 2 (p t ), F 2 ((f) are calculated respectively for 
both quark and antiquark jets. In order to avoid the influence of a non-flat distribution 
of the variables y,p t and <p on the investigation of the dynamical fluctuations, all 
variables are transformed into their corresponding cumulant forms |T0| . 



Then the fitting to the saturated curves for these second order NFM's vs. the par- 
tition number M — 1, 2, . . . , 40 in ID phase space is carried out and the corresponding 
parameters j y , 7 Pt and 7^ are obtained. In order to eliminate the influence of the mo- 
mentum conservation, the first (for F 2 (<f) the first three) point (s) are omitted when 
fitting the data to Eq.(4). 
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Fig. 3 The parameter 7 as function of y C ut- (o) quark jets (b) antiquark jets 

5. Results and discusion 

The variation of the three 7 a (a = y,p t ,f) with y cnt are shown In Fig. 3. From the 
figure some properties of the dynamical fluctuations inside the quark and antiquark 
jets can be extracted. 

The most striking feature in Fig. 3 is the variation of 7^ with the parameter y cnt . 
At small values of y cut , 7^ lies at the bottom. When y cut increases 7^ increases rapidly 
and at a certain point 7^ crosses over j pt , turning from < j pt to 7^ > 7 Pt . A distinct 
point called transition point [[J, where 7^ = 7 Pt , exists in both quark- (Fig.3a) and 
antiquark (Fig.36) jets with only slightly different values of y cnt — 0.0047 for quark 
jets and 0.0050 for antiquark jets. 

At the transition point, the Hurst exponents are 

= 1 + 7pt = 0.70 ± 0.03; 0.74 ± 0.03 

1 + lv 

H w = 1 + 7y = 0.70 ± 0.06; 0.75 ±0.04 

1 + ly 

H ptv = 1 +llp = 1.00 ±0.10; 1.01 ±0.06 

1 + 7 Pt 
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The first values are for quark jets and the second ones for antiquark jets. These 
results indicate that at the transition point the dynamical fluctuations inside quark 
(antiquark) jet are anisotropic in the longitudinal-transverse planes, (y,Pt) and (y, ip), 
and isotropic in the transverse plane (p t , f). This means that the quark and antiquark 
jets are circular in the transverse plane with respect to the dynamical fluctuations at the 
transition point and therefore they are called circular jets 0. For these jets H yPt ^ 1, 
H yip ^ 1 showing that the quark (antiquark) jet system is self-afnne fractal 

In the Durham algorithm used in our calculation the relative transverse momentum 

(7) 



k t is related to y cut as fl2 | 

kt — v^cut ■ \fs 

The value of k t corresponding to the transition point is approximately 6.25 GeV for 
quark jet and 6.46 GeV for antiquark jet. The values of various parameters for both 
quark and antiquark jets at the transition point are listed in Table I 

Table I Parameters 7, Hurst exponent H, cut-parameter y cui and relative transverse 

momentum k t at the transition point 



subsample 


Vent 


ly 


7j»t 


1<p 


Hy pt 




Hp tt p 


h (GeV) 


quark 


0.0047 


1.464 


0.730 


0.735 


0.70 


0.70 


1.00 


6.25 


jet 


±0.0003 


±0.036 


±0.041 


±0.123 


±0.03 


±0.06 


±0.10 


±0.17 


antiquark 


0.0050 


1.409 


0.791 


0.803 


0.74 


0.75 


1.01 


6.46 


jet 


±0.0009 


±0.033 


±0.039 


±0.059 


±0.03 


±0.04 


±0.06 


±0.58 



Having analyzed the dynamical fluctuations inside quark (antiquark) jets, it is in- 
teresting to discuss the following question. In high energy experiments some algorithm, 
e.g. the Durham algorithm, is used to select jets. In doing so, the jet sample obtained 
strongly depends on the cut parameter y cut . When the value of y cut is chosen appro- 
priately, the selected jet will be directly observable in the experiments as a bunch of 
particles. These directly observable jets are refered to as visible jets 0. The question 
is: Is the circular jets determined by the properties of the dynamical fluctuations at 
the transition point consistent with the visible jets in experiments? 

In order to see the relation between the circular jets determined by the transition 
point and the visible jets observed in experiments, let us notice the fact that it was 
in the energy region 17-30 GeV that a third jet (the gluon jet) was historically first 



observed in e + e~ collisions [II]. In Fig.4's are shown the ratios R3 of 3-jet events as 
function of the relative transverse momentum k t at 5 different energies ranging from 
11 to 91.2 GeV. The dashed vertical lines correspond to the transition point fc t =6.25 
GeV and k t = 6.46 GeV, respectively. 
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Fig.4 The ratio of R3 as function of the relative transverse momentum k t 



It can be seen from the figure that for the circular jets defined by the transition 
point k t = 6.25 ~ 6.46 GeV (dashed lines in the figure), 3-jet events become noticeable 
just at y/s 17 GeV. This means that the third (gluon) jet first observed at y/s = 17 
GeV is consistent with being the circular jet. 

4. Conclusion 

In this paper we identify the hadronic quark- and antiquark-jets through matching 
them with the parton-level jets. Then we study the dynamical fluctuations inside the 
quark jets and the antiquark jets respectively using Monte Carlo simulation. 

The second order ID factorial Moments of the three variables y,pt, <f are calculated. 
By fitting the results to the projection formula Eq.(4) the variation of the saturation 
exponents 7 2/ ,7 Pt and 7^ with y cnt is obtained. A transition point exists in both the 
quark and the antiquark jets. At the transition point the dynamical fluctuations inside 
the quark jets and the antiquark jets, being similar to each other, are anisotropic in 
the longitudinal-transverse planes and isotropic in the transverse plane. These jets are 
refered to as circular jets. 

The percentage R 3 of 3-jet events for the jets of different k t have been calculated. 
It turns out that for the circular jets defined by the transition point k t = 6.25 ~ 6.46 
GeV, the percentage R s of 3-jet events becomes noticeable just at y/s m 17 GeV, where 
the 3-jet event was historically first observed. Therefore, the circular jets defined by the 
k t at transition point is consistent to be the visible jets, which are directly observable 
in experiments. 
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